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1. Introduction and Review 

This paper is a continuation of the first two parts of this series ([I], [II]). It relies 
on the formulation of the Kepler conjecture in [F]. The terminology and notation 
of this paper are consistent with these earlier papers, and we refer to results from 
them by prefixing the relevant section numbers with I, II, or F. Around each vertex 
is a modification of the Voronoi cell, called the T^-cell and a collection of quarters 
and quasi-regular tetrahedra. These objects constitute the decomposition star at 
the vertex. A decomposition star may be decomposed into standard clusters. By 
definition, a standard cluster is the part of the given decomposition star that lies 
over a given standard region on the unit sphere. 

A real number, called the score, is attached to each cluster. Each star receives a 
score by summing the scores o~(R) for the clusters R in the star. The scores are 
measured in multiples of a point (pt), where pt ss 0.055. If every star scores at most 
8 pt, then the Kepler conjecture follows. 

The steps of the Kepler conjecture, as outlined in Part I, are 

1. A proof that if all standard regions are triangular, the total score is less than 
8pt 

2. A proof that the standard regions with more than three sides score at most Opt 

3. A proof that if all of the standard regions are triangles or quadrilaterals, then 

the total score is less than 8 pt (excluding the case of pentagonal prisms) 

4. A proof that if some standard region has more than four sides, then the star 

scores less than 8 pt 

5. A proof that pentagonal prisms score less than 8pt 
This paper completes the third of these steps. 

The standard regions of the decomposition stars in the face-centered cubic and the 
hexagonal-close packings are regular triangles and quadrilaterals. These stars score 
exactly 8pt. The local optimality results for hexagonal-close packings and face- 
centered cubic packings have been established in [II]. If the planar map is that of 
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a close packing, there are eight quasi-regular tetrahedra. To score 8 pt, the quasi- 
regular tetrahedra must be regular of edge length 2 [1.9.1]. Such decomposition 
stars are precisely those of the close packings. Conjecturally, all the decomposition 
stars have scores strictly less than 8pt, so that various approximations may be 
introduced to prove the desired bounds. 

The standard regions of pentagonal prisms are triangles and quadrilaterals (10 
triangles and 5 quadrilaterals). The pentagonal prisms are the subject of the fifth 
step of the program. 

Theorem 1. Let D* be a decomposition star whose combinatorial structure is not a 
pentagonal prism. Suppose that each standard region of D* is a triangle or quadri- 
lateral. Then the score of D* is at most 8 pt. Equality is attained exactly when 
the decomposition of the unit sphere into standard regions coincides with that of a 
decomposition star in the face-centered cubic or hexagonal-close packing. 

The proof of Theorem 1 relies on many computer calculations. We make a list 
of combinatorial properties that a decomposition star must have for it to have a 
possibility of scoring more than 8 pt. We then make a computer search to find all 
decompositions of the unit sphere into triangles and quadrilaterals that satisfy all 
the properties on the list. The computer search produces an explicit list containing 
nearly 2000 combinatorial types. 

For each of these combinatorial arrangements of triangles and quadrilaterals, we 
have a nonlinear optimization problem: maximize the score over the space of all 
decomposition stars D* with the given combinatorial arrangement. It is not neces- 
sary to solve this optimization problem. It is sufficient to establish an upper bound 
of 8pt. To do this, we define a linear relaxation of the original problem, that is, 
a linear programming problem whose solution strictly dominates the global maxi- 
mum of the original nonlinear problem. This gives an upper bound on the linear 
problem, which is usually less than 8pt. 

In some cases, this procedure leads to a bound greater than 8 pt, and further analysis 
will be required. 

One advantage of our method of linear relaxation is that the verification of the 
bounds is particularly simple. If the linear relaxation asks to maximize c- x subject 
to the system of linear inequalities Ax < b and x > 0, then duality theory produces 
a vector z with nonnegative entries such that c < zA. To verify the bound of 
8pt, it is enough to check that c < zA, z > 0, and z ■ b < 8pt (because then 
c • x < zAx < z ■ b < 8 pt) . 

Many of the linear inequalities that were used in the linearly relaxed optimization 
problems were obtained as follows. We use numerical methods to find a convex 
polygon containing the set of ordered pairs 

(dih(S),a(S)) 

as S ranges over quasi-regular tetrahedra. The edges of the polygon give linear 
inequalities relating the dihedral angles to the score a. More generally, additional 
inequalities are obtained by considering polygons that contain the ordered pairs 

(dih(S),cr(S*) - Asol(S*)), 
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for appropriate constants A. Similar inequalities are obtained for quadrilateral 
regions. 

We must remember that ultimately we arc dealing with a nonlinear optimization 
problem that is larger, by a considerable order, than what is conventionally thought 
to be solvable by exact methods. The domain of our optimization problem has 
many components, and the dimension varies from component to component. Even 
the magnitude of the problem is poorly understood. The best-known bound on 
the dimension of the components is about 155 dimensions. Components of interest 
frequently have more than 35 dimensions. The decomposition into standard regions 
gives the problem an aspect of separability. Nevertheless, certain complications will 
have to be tolerated. 

The biggest weakness of this method is that the output from the computer search 
for the combinatorial arrangements of triangles and quadrilaterals is not easily 
checked for errors. The algorithm is described in some detail in Section 9, but the 
only assurance that no cases have been skipped comes through a careful reading of 
the computer code. It would be advantageous to have a more transparent proof of 
the results of this section. 

This paper is supplemented by an appendix giving further details of the combina- 
torial arrangements. Additional details about these calculations, including the full 
source code for all of the computer verifications of this paper, can be found in [H2]. 

The outline of this paper was developed at the University of Chicago during the 
summer of 1994. I would like to thank to P. Sally for making computer resources and 
other facilities available to me during my stay in Chicago, and for his encouragement 
with this work. I would also like to give special thanks to S. Ferguson for many 
helpful discussions concerning this topic. His investigations have led to a number 
of improvements in the results presented here. 

2. Geometric Considerations 



2.1. We will call a standard cluster over a quadrilateral region a quad cluster. The 
four vertices of the quad cluster whose projections to the unit sphere mark the 
extreme points of the quadrilateral region will be called the corners of the cluster. 
We call the four angles of the standard region associated with the quad cluster its 
dihedral angles. 

The rules defining the score have undergone a long series of revisions over the last 
several years. The formulation used in this paper is described in [F]. 

Lemma 2.2. A quadrilateral region does not enclose any vertices of height 2.51 or 
less. 

Proof. Let v\ , . . . ,V4 be the corners of the quad cluster, and let v be an enclosed 
vertex of height at most 2.51. We cannot have |i>, — v\ < 2.51 for two different 
vertices Vi, because two such inequalities would partition the region into two sep- 
arate standard regions instead of a single quadrilateral region. We apply 1.4.3 to 
simplify the quad cluster. (Lemma 1.4.3 assumes the existence of another enclosed 
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vertex v' , but it can be omitted both from the statement of the Lemma and from 
the proofs without affecting matters.) Then 1.4.3 allows us to assume 

\vi - v i+ i\ = 2.51, \vi\=2, \v\ = 2.51, 

for i = 1, ... ,4. Reindexing and perturbing w as necessary, we may assume that 
2 < |i>i — v | < 2.51 and |uj — i>| > 2.51, for i = 2,3,4. Moving v, we may assume 
it reaches the minimal distance to two adjacent corners (2 for vi or 2.51 for Vi, 
i > 1). Keeping v fixed at this minimal distance, perturb the quad cluster along 
its remaining degree of freedom until v attains its minimal distance to three of the 
corners. This is a rigid figure. There are four possibilities depending on which 
three corners are chosen. Pick coordinates to show that the distance from v to the 
remaining vertex violates its inequality. □ 

3. Functions Related to the Score 

Set C' 1 := sol(S'(2,2,2,2,2,2)) = 2 arctan(V2/5). The constant C is related to 
the other fundamental constants by the relations pt — 2/( — ir/3 and S oc t — (ir — 
2/0/V8. Rogers's bound is V2/C « 0.7796. 

We consider the functions <J\(R) := <j{R) — A£sol(i?) pt, for A = 0, 1, or 3.2, where 
R is a standard cluster. The constant 3.2 was determined experimentally. We will 
see that <j\ has a simple interpretation. We write t(R) = —<Ti(R). If D* is a 
decomposition star with standard clusters {R}, set t(D*) = J2r t (R)- 

Lemma 3.1. t(R) > 0, for all standard clusters R. 

Proof. If R is not a quasi- regular tetrahedron, or if it is but rad(i?) > 1.41, then 
a(R) < and sol(i?) > 0, so that the result is immediate (see 1.9.17). Assume 
that R is a quasi-regular tetrahedron and rad(i?) < 1.41. The result follows from 
Calculation 10.3.6, which asserts that T(R) < so\(R)(pt. (Equality is attained only 
for the regular tetrahedron of edge 2.) □ 

Lemma 3.2. 

cr(ZT) = Air(pt-T{D*). 

Proof. Let Ri, . . . , Rk be the standard clusters in D* . Then 

a(D*) = J2°(Ri) + (4?r - ^sol(i?;))C^= ^(pt- ^(ifc). 

□ 

Since A-kQ < 22.8, we find as an immediate corollary that if there are standard 
clusters satisfying r(i?i) + • • • + r(Rk) > 14.8 pt, then the score of the star is less 
than 8 pt. 

The function t(R) gives the amount squandered by a particular standard cluster R. 
If nothing is squandered, then r(Ri) = for every standard cluster, and the upper 
bound is AirQpt s=s 22.8 pt. This is Rogers's bound on density. It is the unattainable 
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bound that would be obtained by packing regular tetrahedra around a common 
vertex with no distortion and no gaps. (More precisely, in the terminology of [HI], 
the score s — 4:ir( pt corresponds to the effective density 167r5 oct /(167r — 3s ) = 
•s/2/C ~ 0.7796, which is Rogers's bound.) Every positive lower bound on r(i?j) 
translates into an improvement on Rogers's bound. To say that a decomposition 
star scores at most 8 pt is to say that at least (47r£ — 8)pt w 14.8 pt are squandered. 

4. Some Linear Constraints 

This section gives some linear inequalities between cr(R) — Xso\(R)(pt and dih(i?). 

Proposition 4.1. Let R be a quad cluster. Let cr(R) denote its score, let dih(i?) be 
one of the four dihedral angles of R, and letsol(R) be the solid angle of the standard 
region of R. The following inequalities hold among dih(i?) ; so^S"), and a{R): 

1: a{R) < -5.7906 + 4.56766 dih(i?), 

2: a(R) < -2.0749 + 1.5094 dih(i?), 

3: o-(R) < -0.8341 + 0.5301 dih(i?), 

4: o-(R) < -0.6284 + 0.3878 dih(i?), 

5: o-(R) < 0.4124- 0.1897 dih(i?), 

6: cr(i?) < 1.5707 - 0.5905 dih(i?), 

7: o-(R) < 0.41717 - 0.3 sol(R), 

8: a!(R) < -5.81446 + 4.49461 dih(i?), 

9: ai(R) < -2.955 + 2.1406 dih(i?), 

10: o-i (R) < -0.6438 + 0.316 dih(i?), 

11: ai(R) < -0.1317, 

12: ai(R) < 0.3825- 0.2365 dih(iJ), 

13: o-i (R) < 1.071 - 0.4747 dih(i?), 

14: cr 3 . 2 (i?) < -5.77942 + 4.25863 dih(i?), 

15: cr 3 . 2 (i?) < -4.893 + 3.5294 dih(ii), 

16: (T 3 .2(i?) < -0.4126, 

17: £7 3 .2 (i?) < 0.33 - 0.316 dih(iJ), 

18: a(R) < -0.419351 sol(i?) - 5.350181 + 4.611391 dih(i?), 
19: a{R) < -0.419351 sol(R) - 1.66174 + 1.582508 dih(iJ), 
20: o-(R) < -0.419351 sol(R) + 0.0895 + 0.342747 dih(R), 
21: o-(R) < -0.419351 sol(i?) + 3.36909 - 0.974137 dih(iJ). 

Proposition 4.2. Let R be a quad cluster. Let dihi(i?) and dih 2 (i?) be two adja- 
cent dihedral angles of R. Set d(R) — dihi (R) + dih 2 (R) . The following inequalities 
hold between d(R) and a(R). 

1: a(R) < -9.494 + 3.0508 d(R), 
2: a{R) < -1.0472 + 0.27605 d(R), 
3: o-(R) < 3.5926- 0.844 d(R), 
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1: 1.153 < dih(i?), 
2: dih( J R) < 3.247. 

Proof. Proposition 4.3 follows from interval arithmetic calculations based on the 
methods of [I]. The lower bound of 1.153 holds in fact for the dihedral angles of 
any standard cluster other than quasi-regular tetrahedra. □ 

All of these inequalities have been proved by interval arithmetic methods by com- 
puter. An appendix gives a general description of the cases involved in the verifi- 
cation. Further details are available [H2]. 

5. Types of Vertices 

The combinatorial structure of a decomposition star is conveniently described as 
a planar map. A planar graph is a graph that can be embedded into the plane or 
sphere. A planar map is a planar graph with additional combinatorial structure 
that encodes a particular embedding of the graph [T]. All our planar maps will 
be unoriented: we do not distinguish between a planar map and its reflection. 
Associated with a planar map are faces, (combinatorial) angles between adjacent 
edges, and so forth. Associated with each planar map L is a planar graph G(L), 
obtained by forgetting the additional combinatorial structure. Each planar map 
has a dual L* , obtained by interchanging faces and vertices. The faces of a planar 
map are in natural bijection with the vertices of L* . We say that a face is an n-gon 
if the corresponding vertex in the dual L* has degree n. The boundary of a face is 
an n-circuit in G(L). The edges of the boundary are in natural bijection with the 
edges in L* that are joined to the dual vertex of the face in L*. 

Associated with each decomposition star is a standard decomposition of the unit 
sphere, as described in Part I. We form a planar map L by associating with each 
standard region a face of L and with each edge of a standard region an edge of 
L. This paper is concerned with the special case of the Kepler conjecture in which 
every face of L is a triangle or quadrilateral. 

We say that a vertex v of L has type (p, q) if there are exactly p triangular faces 
and q quadrilateral faces that meet at v. We write (p v ,q v ) for the type of v. The 
type of a vertex of the decomposition star is the type of the corresponding vertex 
in the planar map. 

We use the following strategy in the proof of step 3 of the Kepler conjecture. The 
linear inequalities that were stated in Section 4 will be combined to give a bound 
on the score of the standard clusters around a given vertex of a given type. This 
bound will depend only on the type of the vertex. The bound comes as the solution 
to the linear programming problem of optimizing the sum of scores, subject to the 
linear constraints of Section 4 and to the constraint that the dihedral angles around 
the vertex sum to 2ir. Similarly, we obtain a lower bound on what is squandered 
around each vertex. 

This gives certain obvious constraints on decomposition stars. For example, if 
more than (AttC, — 8) pt are squandered at a vertex of a given type, then that type 
of vertex cannot be part of a decomposition star scoring more than 8 pt. These 
relations between scores and vertex types will allow us to reduce the feasible planar 
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maps to an explicit finite list. For each of the planar maps on this list, we calculate 
a second, more refined linear programming bound on the score. Often, the refined 
linear programming bound is less than 8 pt. 

This section derives the bounds on the scores of the clusters around a given vertex 
as a function of the type of the vertex. Define constants T^p(p, q)/pt by Table 5.1. 
The entries marked with an asterisk will not be needed. 



(5.1) 



t lp(p, o)/pt 


q = 


1 


2 


3 


4 


5 


p = 


* 


* 


15.18 


7.135 


10.6497 


22.27 


l 


* 


* 


6.95 


7.135 


17.62 


32.3 


2 


* 


8.5 


4.756 


12.9814 


* 


* 


3 




3.6426 


8.334 


20.9 






4 


4.1396 


3.7812 


16.11 




* 




5 


0.55 


11.22 




* 




* 


6 


6.339 


* 


* 


* 


* 


* 


7 


14.76 




* 




* 


* 



Proposition 5.2. Let Si,... ,S P and Ri, . . . ,R q be the tetrahedra and quad clus- 
ters around a vertex of type (p, q). Consider the constants of Table 5.1. We have 

V Q 

+ T(ifc) >T LP (p,q), 

Proof. Set 

(d?,t?) = (dih(Si), •>-(#)), {d\,t\) = (dih^), r(fli))- 

Then J2 P r(Si) + J^ 9 T {Ri) is at least the minimum of YT A + subject to 

J2 P d'i + d\ — 2tt and to the system of linear inequalities of Section 10 (Group 3) 
and Proposition 4.1 (obtained by replacing —<j\ and dihedral angles by t\ and d\). 
The constant tlp (p, q) was chosen to be slightly larger than the actual minimum 
of this linear programming problem. 

The entry tlp (5, 0) is based on Lemma 5.3, k = 1. □ 

Lemma 5.3. Let Vi, . . . ,Vk, for some k < 4, be distinct vertices of a decomposi- 
tion star of type (5, 0). Let Si, . . . ,S r be quasi-regular tetrahedra around the edges 
(0, Vi), for i < k. Then 

r 

^2r(Si) > 0.55k pt, 

i=i 

and 

r 

^cr(5i) < rpt- 0A8kpt. 

i=l 
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Proof. We have t(S) > 0, for any quasi-regular tetrahedron S. We refer to the 
edges 7/4, ?/5, ye of a simplex S(yi, . . . ,y&) as its top edges. Set £ = 2.1773. 

We claim (Claim 1) that if Si, . . . , S5 are quasi-regular tetrahedra around an edge 
(0, v) and if Si = S(yi, . . . , y§), where 2/5 > £ is the length of a top edge e on Si 
shared with S2, then J2i T (Si) > 3(0.55) pt. This claim follows from Inequalities 
10.5.1 and 10.5.2 if some other top edge in this group of quasi-regular tetrahedra 
has length greater than £. Assuming all the top edges other than e have length at 
most £, the estimate follows from J2i dih(S,) = 2n and Inequalities 10.5.3, 10.5.4. 

Now let Si, . . . , Ss be the eight quasi-regular tetrahedra around two edges (0, i>i), 
(0,^2) of type (5,0). Let Si and S2 be the simplices along the face (0, vi,V2)- 
Suppose that the top edge (^1,^2) has length at least £. We claim (Claim 2) that 
J2i T (Si) > 4(0.55) If there is a top edge of length at least £ that does not lie 
on Si or S2, then this claim reduces to Inequality 10.5.1 and Claim 1. If any of 
the top edges of Si or S2 other than (vi,V2) has length at least £, then the claim 
follows from Inequalities 10.5.1 and 10.5.2. We assume all top edges other than 
(vi,V2) have length at most £. The claim now follows from Inequalities 10.5.3 and 
10.5.5, since the dihedral angles around each vertex sum to 2n. 

We prove the bounds for r. The proof for a is entirely similar, but uses the constant 
£ = 2.177303 and the Inequalities 10.5.8-10.5.14 rather than 10.5.1-10.5.7. Claims 
analogous to Claims 1 and 2 hold for the a bound by Inequalities 10.5.8-10.5.12. 

Consider r for k = 1. If a top edge has length at least £, this is Inequality 10.5.1. If 
all top edges have length less than £, this is Inequality 10.5.3, since dihedral angles 
sum to 2tt. 

We say that a top edge lies around a vertex v if it is an edge of a quasi-regular 
tetrahedron with vertex v. We do not require v to be the endpoint of the edge. 

Take k — 2. If there is an edge of length at least £ that lies around only one of 
Vi and V2, then Inequality 10.5.1 reduces us to the case k = 1. Any other edge of 
length at least £ is covered by Claim 1 . So we may assume that all top edges have 
length less than £. And then the result follows easily from Inequalities 10.5.3 and 
10.5.6. 

Take k = 3. If there is an edge of length at least £ lying around only one of the Vi, 
then Inequality 10.5.1 reduces us to the case k = 2. If an edge of length at least 
£ lies around exactly two of the v i7 then it is an edge of two of the quasi-regular 
tetrahedra. These quasi-regular tetrahedra give 2(0.55) pi, and the quasi-regular 
tetrahedra around the third vertex Vi give 0.55 pt more. If a top edge of length at 
least £ lies around all three vertices, then one of the endpoints of the edge lies in 
{vi,V2,Vs}, so the result follows from Claim 1. Finally, if all top edges have length 
at most £, we use Inequalities 10.5.3, 10.5.6, 10.5.7. 

Take k = 4. Suppose there is a top edge e of length at least £. If e lies around 
only one of the Vi, we reduce to the case k = 3. If it lies around two of them, 
then the two quasi-regular tetrahedra along this edge give 2(0.55) pt and the quasi- 
regular tetrahedra around the other two vertices Vi give another 2(0.55) pt. If both 
endpoints of e are among the vertices Vi, the result follows from Claim 2. This 
happens in particular if e lies around four vertices. If e lies around only three 



SPHERE PACKINGS III 



9 



vertices, one of its endpoints is one of the vertices Vi, say v\. Assume e is not 
around v 2 . If v 2 is not adjacent to v\, then Claim 1 gives the result. So taking v\ 
adjacent to v 2 , we adapt Claim 1, by using Inequalities 10.5.1-10.5.7, to show that 
the eight quasi-regular tetrahedra around Vi and v 2 give 4(0.55) pt. Finally, if all 
top edges have length at most £, we use Inequalities 10.5.3, 10.5.6, 10.5.7. □ 



Recall that a vertex of a planar map has type (p, q) if it is the vertex of exactly p 
triangles and q quadrilaterals. This section restricts the possible types that appear 
in a decomposition star. 

Let ti denote the constant 0.1317 w 2.37838774 Proposition 4.1.11 asserts that 
every quad cluster R satisfies t(R) > t\. 

Lemma 6.1. The following eight types (p,q) are impossible: (1) p > 8, (2) p > 6 
and q > 1, (3) p > 5 and q > 2, (}) p > 4 and q > 3, (5) p > 2 and q > 4, (6) 
p > and q > 6, (7) p < 3 and q = 0, (8) p < 1 and q = 1. 

Proof. By Proposition 4.1.3 and Calculation 10.1.3, a lower bound on the dihedral 
angle of p simplices and q quadrilaterals is 0.8638p + 1.153g. If the type exists, 
this constant must be at most 2n. One readily verifies in Cases 1-6 that 0.8638p + 
1.153g > 2ir. By Proposition 4.3 and Calculation 10.1.2, an upper bound on the 
dihedral angle of p triangles and q quadrilaterals is 1.874445p+ 3.247q. In Cases 7 
and 8 this constant is less than 2n. □ 

Lemma 6.2. // the type of any vertex of a decomposition star is one of (4,2), 
(3,3), (1,4), (1,5), (0,5), (0,2), (7,0), then the decomposition star scores less than 
8pt. 

Proof. According to Table 5.1, we have tlp(p, q) > (47r£ — 8) pt, for (p, q) = (4, 2), 
(3, 3), (1, 4), (1, 5), (0, 5), or (0, 2). By Lemma 3.2, the result follows in these cases. 
Now suppose that one of the vertices has type (7,0). By the results of Part I, 
which treats the case in which all standard regions are triangles, we may assume 
that the star has at least one quadrilateral. We then have t(D*) > tlp(7, 0) +t/ L > 
(4ttC - 8) pt. The result follows. □ 

In summary of the preceding two lemmas, we find that we may restrict our attention 
to the following types of vertices. 



6. Limitations on Types 



(6,0) 
(5,0) 
(4,0) 



(5,1) 
(4,1) 
(3,1) 
(2,1) 



(3,2) 

(2,2) (2,3) 
(1,2) (1,3) 



(0,3) (0,4) 
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7. Properties of Planar Maps 

Proposition 7.1. Suppose that cr(D*) > 8pt. The planar map L of D* has the 

following properties (without loss of generality): 

1. The graph G(L) has no loops or multiple joins. 

2. Each face of L is a triangle or quadrilateral. 

3. L has at least eight triangular faces. 

4. L has at most six quadrilateral faces, and at least one. 

5. Each vertex has one of the following types: (6,0), (5,0), (4,0), (5,1), (4,1), 
(3,1), (2,1), (3,2), (2,2), (1,2), (2,3), (1,3), (0,3), and (0,4). 

6. If C is a 3-circuit in G(L), then it bounds a triangular face. 

7. If C is a 4-circuit in G(L), then one of the following is true: 

(a) C bounds some quadrilateral region, 

(b) C bounds a pair of adjacent triangles, 

(c) C encloses one vertex, and it has type (4,0) or (2, 1). 

8. t^q + J2v£v( t Lp(P v >1 v ) ~ t*Qv) < (4ttC — 8) pt, for any collection V of vertices 
in L such that no two vertices of V lie on a common face. 

Proof. A loop would give a closed geodesic on the unit sphere of length less than w. 
A multiple join would give nonantipodal conjugate points on the sphere. Property 
2 is the restriction of step 3 of the Kepler conjecture. Property 3 follows from 
1.9.1 and Part II. Property 4 follows from lt^ > (4-7r£ — 8) pt. If there are no 
quadrilaterals, the problem has been solved in Part I. The restrictions on types 
were obtained in Section 6. Property 6 is established in Part I. A 4-circuit encloses 
at most one vertex by 1.4.2. If it encloses none, it gives a quad cluster or two 
tetrahedra. Otherwise, it encloses a vertex of type (4,0) or type (2, 1). Property 8 
comes from Section 5. □ 

8. Combinatorics 

In Steps III and IV of the Kepler conjecture, we need to generate all planar maps 
satisfying various lists of conditions. Here we describe a computer algorithm, which 
has been implemented in Java. We describe the algorithm in a way that it can be 
used for Step IV as well. 

We assume that the planar maps satisfy the following conditions. 

1. There are no loops or multiple joins. 

2. Each face is a polygon. 

3. The graph has between 3 and N vertices, for some explicit N. 

4. The degree at each vertex is at most 6. 

5. If C is a 3-circuit in G(L), then C bounds a triangular face. 
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6. If C is a 4-circuit in G(L), then one side of C contains at most 1 vertex. 

7. There is a constant T > and constants t n > and t(p, q, r) > 0, such that 



where the first sum runs over finished faces, and the second sum runs over any 
separated set of vertices. (The term finished will be described below. For planar 
maps produced as output of the algorithm, every face will be finished. But at 
intermediate stages of the algorithm some will be unfinished.) Also, a separated set 
of vertices means a collection of vertices with the properties that no two lie on a 
common face, and that every face at each of the vertices is finished. Here p v is the 
number of triangles, q v is the number of quadrilaterals, r v is the number of other 
regions at v, and n = n v = p v + q v + r v . Also, n is the number of sides of the 
finished face. 

There are additional properties that it might be helpful to impose, but the ones 
stated are sufficient for the description of the algorithm. In our context, r = r v = 0, 
t 3 = 0, T = (4tt< - 8) pt, U = 0.1317 and t(p, q, 0) = tl P (p, q) - t 4 q. 

We produce all maps satisfying the conditions 1-7, by extending maps already 
satisfying these properties. We assume that we have a stack of maps {Lj : j e J} 
satisfying the conditions 1-7 such that each face of each of these maps is labeled 
finished or unfinished. Each planar map in the stack will have at least one unfinished 
face. In one iteration of the algorithm, we pop one of the maps L from the stack, 
modify it in various ways to produce new maps, output any of the new maps that 
are finished (meaning all faces are finished), and push the remaining ones back on 
the stack. When all maps have been popped from the stack, we are guaranteed to 
have produced all finished maps satisfying properties (1-7). Here are the details of 
the algorithm. 

1. Let L be a planar map that has been popped from the stack. Fix any unfinished 
face F of L and any edge e of F. Label the vertices of F consecutively 1, . . . ,£, 
with 1 and t the endpoints of e. For each m = 3, . . . , A, let A m be the set of all 
m-tuples (cii) € Z m satisfying (1 = a\ < < ■ ■ ■ < a m = £), with a m _i ^ a m . 

2. For each a G A m , wc draw a new m-gon F a along the edge e of F as follows. We 
construct the vertices v(l), . . . , v(m) inductively If i = 1 or a, ^ ai-i, then we set 
v(i) = vertex dj of F. But if a, = aj_i, we add a new vertex v(i) to the planar 
map. The new face is to be drawn along the edge e over the face F. (For example, 
when I = 5, the faces F a corresponding to a = (1,1,3,4,4,5) and a = (1,1,1,5) 
are shown in Diagram 8.1.) As we run over all m and a e A m , wc run over all 
possibilities for the finished face along the edge e inside F. 




i 



v 
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Diagram 8.1 



3 



3 



4 




2 



3. The face F a is to be marked as finished. By drawing F a , F is broken into a 
number of smaller polygons. (In Diagram 8.1, F is replaced respectively by three 
and two polygons.) Each of these smaller polygons other than F a is taken to be 
unfinished, except triangles, which are always taken to be finished. 

4. Various planar maps extending L are obtained by this process. Those that do 
not satisfy the conditions (1-7) are discarded. Those that have no unfinished faces 
are output. The remaining ones are pushed back onto the stack. If the stack is 
empty, the program terminates. Otherwise, we pop a planar map from the stack 
and return to the first step. It is condition 7 that forces the algorithm to terminate. 

To begin the algorithm, we need an initial stack of planar maps. The planar maps 
in the initial stack will be called seeds. To produce a list of seeds, it is enough to give 
any list that is guaranteed to produce all possibilities by the algorithm described 
in 1, . . . ,4. For example, for Part III we want all configurations with at least one 
quadrilateral and nothing but triangles and quadrilaterals. We could let our initial 
stack consist of a single planar map L, where the graph G(L) is a single 4-cycle. 
Also, L is to have two faces, a finished quadrilateral and an unfinished complement. 
Then by iterating through the steps 1, . . . , 4, we generate all possible extensions of 
a quadrilateral to a planar map satisfying 1, . . . ,7. 

Although this seed would work, in order to improve the performance of the algo- 
rithm, in the implementation used for this paper, we used a more detailed list of 
seeds, based on the classification of types (p, q) in Section 6. 

This algorithm produced 1762 cases, even when the additional properties listed in 
Proposition 7.1 were used. To be exact, a few of the maps may be superfluous, 
because there was no need to discard every last map that we were allowed to. The 
important point is that an explicit finite list was obtained. Because of the number 
of possibilities involved we have not listed them here. The Java source code and 
pictures of the maps are available at [H2] . 



9. Linear Programming Bounds 



For each of the planar maps produced in Section 8, we define a linear programming 
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problem whose solution dominates the score of the decomposition stars associated 
with the planar map. A description of the linear programs is presented in this 
section. 

Theorem 9.1. Let L be any planar map obtained in Section 8. One of the following 
holds. (1) L is the planar map of the pentagonal prism, hexagonal-close packing, 
or face- centered cubic. (2) Every decomposition star with planar map L scores less 
than 8pt. (3) L is one of the 18 cases presented in Appendix I. 

The variables of the linear programming problem are the dihedral angles, the scores 
of each of the standard clusters, and their edge lengths. 

We subject these variables to a system of linear inequalities. First of all, the dihedral 
angles around each vertex sum to 2n. The dihedral angles, solid angles, and score 
are related by the linear inequalities of Groups 1, 2, 3, and 4 in Section 10. These 
include Propositions 4.1 and 4.2. The solid-angle variables are linear functions of 
dihedral angles. The score of a decomposition star is 



Forgetting the origin of the scores, solid angles, and dihedral angles as nonlinear 
functions of the standard clusters and treating them as formal variables subject 
only to the given linear inequalities, we obtain a linear programming bound on the 
score. 

Floating-point arithmetic was used freely in obtaining these bounds. The linear 
programming package CPLEX was used (see www.cplex.com). However, the results, 
once obtained, could be checked rigorously as follows. (We did not actually do 
this because the precision never seemed to be an issue, but this is how it can 
easily be done. 1 ) For each quasi- regular tetrahedron Si we have a nonnegativc 
variable Xi = pt— cr(Si). For each quad cluster i?^, we have a nonnegative variable 
Xk = —a(Rk). A bound on the score is ppt— ^2 ieI Xi, where p is the number of 
triangular standard regions, and / indexes the faces of the planar map. We give 
error bounds for a linear program involving scores and dihedral angles. Similar 
estimates can be made if there are edges representing edge lengths. Let the dihedral 
angles be Xj, for j in some indexing set J. Write the linear constraints as Ax < b. 
We wish to maximize c • x subject to these constraints, where Cj = — 1, for i e J, 
and Cj = 0, for j G J. Let z be an approximate solution to the inequalities zA > c 
and z > obtained by numerical methods. Replacing the negative entries of z by 
we may assume that z > and that zA{ > Cj — e, for i G / U J, and some small 
error e. If we obtain the numerical bound ppt + z ■ b < 7.9999 pt, and if e < 1CP 8 , 
then the score is less than 8pt. In fact, note that 



is at least Cj for i e / (since Q = —1), and that it is greater than Cj — e/(l + e), for 
i G J (since Ci = 0). Thus, if N < 60 is the number of vertices, and p < 2(N — 2) < 

lr Thc output from each linear program in this paper has been double checked with interval arith- 
metic. Predictably, the error bounds presented here were satisfactory. 1/2002 



<r(5i) + • • • + o-(S p ) + <7(i?i) + • • • + <j(R q ). 
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116 is the number of triangular faces, 



a(D*) < ppt + c • x < ppt + 



1 + e 



z 



) 




z ■ o e „ „ T 
< ppt+ h 2?riV 



1+e 1+e 



< [ppt+ z-b + e(ppt+2nN)]/(l + e) 

< [7.9999 pt+ 10" 8 (116pt+ 500)] /(l + 10~ 8 ) < 8 pi 



In practice, we used 0.4429 < 0.79984 pt as our cutoff, and N < 14 in the interesting 
cases, so much tighter error estimates are possible. 



In each of these calculations, when the cluster is a quasi-regular tetrahedron 5, we 
set a — cr(S), dih = dih(S'), and so forth, Let a\ = a — A£ptsol, for A = 1, 3.2. We 
make similar abbreviations for quad clusters. The inequalities in Group 1 follow 
from results appearing elsewhere. These inequalities have been verified by interval 
arithmetic in [H2]. 

Group 1. Calculations that have been verified elsewhere. 

Quasi-regular tetrahedra: 

1. cr < pt (1.9.1), 

2. dih < 1.874445 (1.8.3.2), 

3. dih > 0.8638 (1.9.3), 

4. cr < -0.37642101 sol +0.287389 (1.9.8), 

5. cr < 0.446634 sol -0.190249 (1.9.9), 

6. cr < -0.419351 sol +0.2856354 + 0.001 (1.9.10,1.9.11,1.9.12,1.9.18). 

Quad clusters: 

7. cr < (II). 



Group 2. Inequalities for quasi-regular tetrahedra depending on edge lengths. 

1. sol > 0.551285 + 0.199235(?/4 + 2/5 + Ve - 6) - 0.377076(?/i + y 2 + y 3 - 6), 

2. sol < 0.551286 + 0.320937(?/4 + 2/5 + Ve - 6) - 0.152679(yi + y 2 + y 3 - 6), 

3. dih > 1.23095-0.359894(t/ 2 + 2/3 + y5 + 2/6-8) + 0.003(t/i-2) + 0.685(y 4 -2), 

4. dih < 1.23096-0.153598(2/2+t/3+y5+2/6-8)+0.498(yi-2)+0.76446(y 4 -2), 

5. cr < 0.0553737 - 0.10857(yi + • • • + y 6 - 12), 

6. cr + 0.419351 sol < 0.28665 - 0.2 (j/i + y 2 + 2/3 - 6), 

7. u\ < IQ- 6 - 0.129119(2/4 + 2/5 + ye - 6) - 0.0845696(2/1 + 2/2 + 2/3 - 6). 



Group 3. General inequalities for quad clusters and quasi-regular tetrahedra. 

Quasi-regular tetrahedra: 

1. cr < 0.37898 dih -0.4111, 

2. cr < -0.142 dih +0.23021, 

3. cr < -0.3302 dih +0.5353, 

4. o-i < 0.3897 dih -0.4666, 



10. Calculations 
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5. 01 < 0.2993 dih -0.3683, 

6. (7i < 0, 

7. cri < -0.1689 dih +0.208, 

8. cri < -0.2529 dih +0.3442, 

9. cr 32 < 0.4233 dih -0.5974, 

10. cr 32 < 0.1083 dih -0.255, 

11. cr 3 2 < -0.0953 dih -0.0045, 

12. cr 3 2 < -0.1966 dih +0.1369, 

13. cr < -0.419351 sol+0.796456dih-0.5786316, 

14. cr < -0.419351 sol +0.0610397 dih+0.211419, 

15. cr < -0.419351 sol-0.0162028dih+0.308526, 

16. cr < -0.419351 sol -0.0499559 dih +0.35641, 

17. cr < -0.419351 sol-0.64713719dih+l. 3225. 
Quad clusters: Propositions 4.1 and 4.2. 



Group 4. Miscellaneous inequalities. 

Quasi-regular tetrahedra: 

1. The quasi-regular tetrahedra at a vertex of type (4, 0) 
score at most 0.33 pt (1.5.2). 

2. The sum of the dihedral angles around a vertex is 2n. 

3. The five quasi-regular tetrahedra Si at a vertex of type (5, 0) 
satisfy 

^2<r{Si) < ^(-0.419351 sol(5i) +0.2856354) 

(1.5.1.1). 

Flat quarters: 

4. -0.398(y 2 + Vs + 2/5 + 2/6 ) + 0.32572/1 - dihi < -4.14938, 
if y 4 > 2.51. 

5. Proposition 4.3, Lemma 5.3. 

6. Inequalities of Appendix 1: A.2.1-11, A.3.1-11, A.4.1-4, A.6.1-9, 
A6.1'-8', A.8.1-3. 



Group 5. Inequalities used by Lemma 5.3. 

quasi-regular tetrahedra: Let £ = 2.1773, m = 0.2384. 

1- If 2/4 > £, then r > 0.55 pt, 

2- If 2/4, 2/5 > £, then r > 2(0.55) pt, 

3- If 2/4 < £, then t > -0.29349 + m dih, 

4. If 2/4, 2/6 < £, 2/5 > £, then r > -0.26303 + m dih, 

5- If 2/6 > 2/4, 2/5 < C, then r > —0.5565 + m(dihi + dih 2 ), 

6. If y A , 2/5, 2/6 < £, then r > -2(0.29349) + m(dihi + dih 2 ), 

7. If ?/4, 2/5, 2/6 < C, then t > -3(0.29349) + m(dihi + dih 2 + dih 3 ). 
Now set = 2.177303, m = 0.207045. 

8- If 2/4 > £, then cr < (1 - 0.48) pt, 

9- If 2/4, 2/5 > then cr < (1 - 2(0.48)) pt, 

10. If ?/4 < £, then ct < 0.31023815 - mdih, 

11. If 2/4,2/6 < 2/5 > £, then cr < 0.28365 - mdih, 

12. If 2/6 > 2/4, 2/5 < I, then a < 0.53852 - m(dihj + dih 2 ), 

13. If 2/4, 2/5, 2/6 < £, then cr < -pt + 2(0.31023815) - m(dihi + dih 2 ), 

14. If 2/4,2/5,2/6 < h then cr < -2 pt + 3(0.31023815) - m(dihi + dih 2 + dih 3 ). 
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Appendix 1. Some Final Cases 

The graphs are numbered in the archive from to 1761. There are three cases that 
are treated elsewhere: PM(4, 1679) is the face-centered cubic, PM(4, 1672) is the 
pentahedral prism, and PM(4, 1640) is the hexagonal close packing. 

The body of this paper eliminates all but a couple of dozen planar maps. (The 
explicit list appears in the archive [H2].) This appendix indicates how to eliminate 
the final cases. The archive contains a few graphs that are isomorphic to each other. 
The following discussion assumes that these duplicates have been eliminated. To 
exploit the nonlincarities of the problem, we use a branch-and-bound method and 
divide the domains of the optimization problem into several thousand smaller sets. 
A linear programming bound of 8 pt is obtained in each case. This appendix lists 
all of the inequalities that have been used, and gives a description of the cases. We 
refer the reader to [H2] for details about computer implementation of the linear 
programs. 

A.l. Types. Each quad cluster with corners (vi, t>2, V3, V4) is one of four types [F]. 
Although it was advantageous to group these cases together to simplify the combi- 
natorics, it is now better to separate these cases and to develop linear inequalities 
for each case. 

1. Two flat quarters with diagonal (i>i, 1*3). The score of each quarter is compression 
or the analytic Voronoi function. 

2. Two flat quarters with diagonal (v2,Vi). The score of each quarter is compression 
or the analytic Voronoi function. 

3. Four upright quarters forming an octahedron. The score of each upright quarter 
is compression or the averaged analytic Voronoi function octavor(Q) = (vor(Q) + 
vor(Q))/2. In other words, <r(Q) = (fi(Q) + n{Q))/2. We are using notation from 
[F]- 

4. One of various mixed quad clusters. The score is at most vor , the truncated 
Voronoi function at radius t = 1.255. 

A. 2. Flat quarters. For flat quarters, we have the following inequalities that 
were established by interval arithmetic. The edge 2/4 is taken to be the diagonal of 
the flat quarter. Here a = fi. 

1. - dih 2 +0.35y 2 - 0.15yi - 0.15j/ 3 + 0.7022y 5 - 0.17y 4 > -0.0123, 

2. - dih 3 +0.35?/3 - 0.15j/i - 0.15y 2 + 0.7022y 6 - 0.17j/ 4 > -0.0123, 

3. dih 2 -0.13|/2 + 0.631J/1 + 0.31j/ 3 - 0.58y 5 + 0.4 1 32/4 + 0.025y 6 > 2.63363, 

4. dih 3 -0.132/3 + 0.6312/1 + 0.31y 2 - 0-58y 6 + 0.4 1 32/4 + 0.025y 5 > 2.63363, 

5. - dihi +0.7142/1 - 0.2212/2 - 0.2212/3 + 0.92j/4 - 0.2212/5 - 0.22l2/ 6 > 0.3482, 

6. dihi -0.3152/1+0.39722/2+0.39722/3-0.7152/4+0.39722/5+0.39722/6 > 2.37095, 

7. -sol -0.1872/1-0. 1872/2 - 0.1872/3 + 0.11852/4 + 0.4792/5 + 0.4792/6 > 0.437235, 

8. sol +0.488j/i + 0.4882/2 + 0+882/3 - 0.334y 5 - 0.3My 6 > 2.244, 

9. -a - 0.1592/1 - O.O8I2/2 - 0.081j/ 3 - 0.133j/ 5 - 0.133y 6 > -1.17401, 

10. cr < -0.419351 sol +0.1448 + 0M36(y 5 + y 6 - 4) + 0.079431 dih, 

11. <t< HT 6 - 0.197(2/4 + 2/5 + 2/6- 2 v/2- 4). 

A. 3. Upright quarters. The following inequalities for upright quarters have 
been established by interval arithmetic. The first edge is taken to be the upright 
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diagonal. 

1. dihi -0.636yi + 0.462y 2 + 0.462y 3 - 0.82y 4 + 0.462y 5 + 0.462y 6 > 1.82419, 

2. - dihi +0.55yi - 0.214j/ 2 - 0.214y 3 + 1.24y 4 - 0.214j/ 5 - 0.214j/ 6 > 0.75281, 

3. dih 2 +0.4yi - 0.15y 2 + 0.09j/ 3 + 0.631y 4 - 0.57y 5 + 0.23y 6 > 2.5481, 

4. - dih 2 -0.454yi + 0.34j/ 2 + 0.154y 3 - 0.346y 4 + 0.805y 5 > -0.3429, 

5. dih 3 +0.4yi - 0.15y 3 + 0.09y 2 + 0.631?/ 4 - 0.57y 6 + 0.23y 5 > 2.5481, 

6. - dih 3 -0.454yi + 0.34j/ 3 + 0.154y 2 - 0.346y 4 + 0.805y 6 > -0.3429, 

7. sol +0.065y 2 + 0.065y 3 + 0.061y 4 - 0.115y 5 - 0.115y 6 > 0.2618, 

8. - sol -0.293yi - 0.03y 2 - 0.03?y 3 + 0.12y 4 + 0.325y 5 + 0.325y 6 > 0.2514, 

9. -a - 0.054y 2 - 0.054y 3 - 0.083y 4 - 0.054j/ 5 - 0.054y 6 > -0.59834, 

10. a < -0.419351 sol +0.079431 dih2 + 0.06904 - 0.0846(j/i - 2.8), 

11- If 2/2, 2/3 < 2.13, then cr < 0.07(j/i - 2.51) - 0.133(y 2 + y 3 + y 5 + y e - 8) - 
0.135(y 4 - 2). 

A. 4. Truncated quad clusters. Let <j>{h,t) = (4 - 25 oc tht(h + t))/3. Set t = 
1.255 and (j>o = 0(£o,^o)- In the truncated case vor , [F] gives 

vor = (j>o sol + ^ MVi/" 2 ) di^ -4<S oct ^ q u0 (- R )> 

R 

with c/) = cj)(to,t ), and 

A(h) = (1 - h/t )(<t>(h, t ) - <p(t , t )). 

Let R be the Rogers simplex i?(z/i/2, 77(2/1, 2/2, ye), to)- The function quo(i?) is de- 
fined in [F.3.3]. We have quo(i?) > 0. Let votq denote the truncated Voronoi score 
of half the quad cluster, divided into two simplices along a diagonal, obtained by 
applying the formula for vor to the simplex. The following inequalities hold by 
interval arithmetic: 

1. dih -0.372j/i + 0.465y 2 + 0.465y 3 + 0.4657/5 + 0.465t/ 6 > 4.885, 

2. - vorf -0.062/ 2 - 0.06j/3 - 0.1852/ 5 - 0.185y 6 > -0.9978, 
provided dih < 2.12, and 2/1,2/2,2/3 < 2.26, 

3. - vorfi +0.419351 so\ A < 0.3072, 
provided dih < 2.12, and 2/1,2/2,2/3 < 2.26, 

4. quo +0.007582/1 + 0.0115t/ 2 + 0.01157/6 > 0.06333. 

Also, A > 0, A' < 0, and A" > 0, for h e [l,t ]- If dih e [dih min , dih max ], and 
h e [1,/imax], for some constants dih m i n , dih max , and 1 < /i max < to, then setting 
A = (A(h max ) — A(l))/(/i max — 1), we obtain the additional elementary inequalities 
for Ad := A(h) dih. 

5. Ad - A(l) dih < A(/i- l)dih min , 

6. Ad < (A(l) + X(h - 1)) dih max . 

We use linear programming methods to determine bounds ft. max , dih m i n , dih max . If 
Ax < b is the system of inequalities used in the linear programs in the main body 
of the paper, then we obtain an upper bound on a variable y by solving the linear 
program max y subject to the constraints Ax < 6, and the constraint that the sum 
of the variables tr (that is, the linear variables corresponding to the score) is at least 
8pt. 
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A. 5. Linear programs. 

Consider one of the remaining cases PM = PM(4, n). Suppose PM has r quadrilat- 
eral faces. We run 4 r linear programs, depending on which type 1-4 of quad cluster 
from A.l each quadrilateral face represents. In each case, we add the additional 
linear inequalities from A. 2, A. 3, or A. 4 as appropriate. Note that a few of these 
inequalities are only conditionally true, so that the inequality can be used only if it 
is known that the condition holds. All of the planar maps have bounds under 8 pt 
by this method, except for PM(4, 71), PM(4, 118), PM(4, 126), and PM(4, 178). 

A. 6. Quasi-regular tetrahedra. 

The first nine inequalities assume that 2/4 + 2/5 + ya < 6.25. Of these, the last two 
are established only under the additional assumptions 2/1,2/2,2/3 < 2.13. The next 
eight inequalities assume that 2/4 + 2/5 + 2/6 > 6.25. Of these, the last three assume 
that 2/1,2/2,2/3 < 2.13. 

1. sol +0.377076j/i + 0.37 7 762/2 + 0.37 7 762/ 3 - 0.221y 4 - O.22I2/5 - 0.221j/ 6 > 
1.487741, 

2. 0.221j/4 + O.22I2/5 + 0.2212/6 - sol > 0.76822, 

3. dih +0.34j/2 + 0.342/3 - 0.689j/ 4 + 0.27y 5 + 0.27y 6 > 2.29295, 

4. - dih +0.4982/1 + 0.7312/4 - 0.2122/5 - 0.212y 6 > 0.37884, 

5. -a - 0.109t/i - 0.1092/2 - 0.109y 3 - 0.14 1 352/4 - 0.14 1 35y 5 - 0.14135y 6 > 
-1.5574737, 

6. -a - 0.419351 sol-0.2j/i - O.22/2 - 0.2j/ 3 - 0.0482/ 4 - 0.0482/ 5 - 0.0482/ 6 > 
-1.77465, 

7. t - 0.08456962/1 - 0.08456962/2 - 0.08456962/ 3 - 0. 163y 4 - 0.163y 5 - 0.163y 6 > 
-1.48542, 

8. dih +0.272/2 + 0.272/3 - 0.6892/4 + 0.272/5 + 0.272/6 > 2.01295, if 2/1, 2/2, 2/3 < 2.13. 

9. -(7-0.141352/1-0.141352/2-0.141352/3-0.141352/4-0.141352/5-0.141352/6 > 
-1.7515737, if 2/1,2/2,2/3 < 2.13. 



!'■ 2/4 + 2/5 + 2/6 > 6.25, 

2'. sol +0.3782/1 + 0.3782/2 + 0.3782/3-0.17812/4-0. 17812/5-0.17812/6 > 1.761445, 
3'. -sol-0.171j/i - 0.1712/2 - 0.1712/3 + 0.3405j/4 + 0.34052/5 + 0.34052/ 6 > 
0.489145, 

4'. -a - O.12O82/1 - 0.12082/2 - O.I2O82/3 - 0.0781j/4 - 0.0781j/5 - 0.0781j/ 6 > 
-1.2436, 

5'. -a- 0.419351 sol -0.2yi - 0.22/ 2 - O.22/3 + O.OIO62/4 + O.OIO62/5 + O.OIO62/6 > 
-1.40816, 

6'. sol+0.3562/i + 0.3562/2 + 0.3562/3-0.178l2/4-0. 17812/5-0.17812/6 > 1.629445, 
if 2/1,2/2,2/3 < 2.13, 

7'. -sol -0.254J/1 - 0.2542/2 - 0.2 54j/ 3 + 0.34052/4 + 0M05y 5 + 0.34 052/ 6 > 
-0.008855, if 2/1,2/2,2/3 < 2.13, 

8'. -cr-0.1672/1-0. 1672/2 - 0. 1672/3 - 0.078l2/4-0.078l2/5-0.078l2/6 > -1.51017, 
if 2/1,2/2,2/3 < 2.13. 



A. 7. Branch and Bound methods. 

For each case that failed the tests of A. 5, we use a branch-and-bound method as 
follows. We pick 10 quasi-regular tetrahedra in the configuration. We divide the 
domain into 2 10 cases by imposing the constraint 2/4+2/5+2/6 < 6.25 or 2/4 + 2/5+2/6 > 
6.25 at each quasi-regular tetrahedron. Depending on which constraint is picked, we 
add (to the inequalities already present) the first group 1-7 or the second group l'-4' 
of inequalities. Whenever we write of the branch-and-bound inequalities associated 
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with various faces, we mean these two groups of inequalities. When r faces are used, 
there are 2 r linear programs. Before running these we determine the maximum 
dihedral angles and edge lengths by separate linear programs. If the conditions 
of A. 4. 2, A. 4. 3 are met, we add these equations. Most of the 2 10 linear programs 
give bounds under 8 pt. The only planar maps that still give bounds over 8 pt are 
PM(4, 126) and PM(4, 178). 

A. 8. Final cases. 

The case PM(4, 178) is not difficult. The maximum height of a vertex is less than 
2.13, by linear programming bounds. Add the inequalities from A. 3 and A. 6 that 
assume this condition. With these additional inequalities, the linear programming 
bound is less than 8 pt. 

In the case PM(4, 126), again we find that the vertices have heights less than 2.13, 
and we add the relevant inequalities. The score is less than 8pt unless the quad 
clusters fall into the pure Voronoi case or mixed cases. (That is, there are no flat 
quarters and no octahedra.) In these remaining cases the diagonals of the quad 
clusters are at least 2y/2. 

Fix a vertex of type (4, 1) in the planar map. There are two such vertices to choose 
from, but there is an isomorphism carrying one to the other. Let vi,V2,V3,V4 be 
the corners of one of the quad clusters, with v\ the chosen vertex of type (4, 1). The 
linear programming upper bound on the dihedral angle of the quad cluster along 
(0,«i) is less than 1.694. The linear programming upper bound on |t> 2 | + \v\ — v 2 \ + 
I vi — V4 1 + |^4 1 is less than 8.709. The following inequality, established by interval 
computations, implies that |i>2 — Vi\ <G [2y/2, 2.93]. 

1. dih > 1.694, if 2/1,2/2, 2/3 < 2.13, 2/2 + 2/3 + 2/5+2/6 < 8.709, and y 4 > 2.93. 

This allows us to add the following two interval computations to the set of linear 
programming inequalities. Both hold for simpliccs satisfying 2/1,2/2,2/3 < 2.13, 2/4 € 
[2\/2, 2.93], 2/5,2/6 <= [2,2.51]. The quad cluster is broken in two simpliccs, with 
these inequalities holding on each simplex. 

2. dih 2 +0.592/1 + O.I2/2 + O.I2/3 + 0.552/4 - O.62/5 - 0.12j/ 6 > 2.6506. 

3. dih 2 +0.352/1 - 0.24?/2 + 0.05y 3 + 0.35y 4 - 0.72y 5 - 0.18j/ 6 < 0.47. 

With these additional inequalities, the linear programming bound drops below 8 pt. 
We conclude that all 18 planar maps that have scores below 8 pt. This concludes 
the proof of the main theorem of the paper and the third step in the proof of the 
Kepler conjecture. 
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Appendix 2. Interval Verifications 

We make a few remarks in this appendix on the verification of the inequalities of 
Proposition 4.1 and 4.2. The basic method in proving an inequality f(x) < for 
x e C, is to use computer-based interval arithmetic to obtain rigorous upper bounds 
on the second derivatives: \fij(x)\ < a,ij, for x e C. These bounds lead immediately 
to upper bounds on f(x) through a Taylor approximation with explicit bounds on 
the error. We divide the domain C as necessary until the Taylor approximation on 
each piece is less that the desired bound. 

Some of the inequalities involve as many as 12 variables, such as the octahedral cases 
of Proposition 4.2. These are not directly accessible by computer. We describe some 
reductions we have used, based on linear programming. We start by applying the 
dimension reduction techniques described in 1.8.7. We have used these whenever 
possible. 

We will describe Proposition 4.2 because in various respects these inequalities have 
been the most difficult to prove, although the verifications of Propositions 4.1 and 
4.3 are quite similar. If there is a diagonal of length < 2\/2, we have two flat 
quarters Si and S 2 . The score breaks up into a = o{S\) + o(Si). The simpliccs 
Si and Si share a three-dimensional face. The inequality we wish to prove has the 
form 

a < a(dih(Si) + dih 2 (Si) + dih 2 (S 2 )) + b. 
We break the shared face into smaller domains on which we have 

<j(Si) < a(dih(Si) + dih 2 (Si)) + h, 
<j(S 2 ) < adih 2 (Si) + 6 2 , 

for some 6i,6 2 satisfying b\ + 6 2 < b. These inequalities are six-dimensional verifi- 
cations. 

If the quad cluster is an octahedron with upright diagonal, there are four upright 
quarters Si, . . . , S4. We consider inequalities of the form 

3 

(X t ) a(S l )<Y,^yj(S l ) + a 7 (dih 1 (S i )-7r/2) + J2^dih J (S t ) + b\ 

J"#4 j=2 

If Y^i=\ a% j — 0' J 7^ 4, and < b, then for appropriate e* C {0,1}, these 

inequalities imply the full inequality for octahedral quad clusters. 

By linear programming techniques, we were able to divide the domain of all octa- 
hedra into about 1200 pieces and find inequalities of this form on each piece, giving 
an explicit list of inequalities that imply Proposition 4.2. The inequalities involve 
six variables and were verified by interval arithmetic. 
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To find the optimal coefficients a* by linear programming we pose the linear problem 

maxi 
such that 

X,, i = 1,2,3,4, {S 1 ,S 2 ,S 3 ,S i }eC 

i 
i 

where {Si, S2, S3, S4} runs over all octahedra in a given domain C . The nonlinear 
inequalities Xi are to be regarded as linear conditions on the coefficients a*-, b l , etc. 
The nonlinear functions cr(Si), dih(Si), Uj(Si) are to be regarded as the coefficients 
of the variables a*-, . . . in a system of linear inequalities. There are infinitely many 
constraints, because the set C of octahedra is infinite. In practice we approximate 
C by a large finite set. If the maximum of t is positive, then we have a collection 
of inequalities in small dimensions that imply the inequality for octahedral quad 
clusters. Otherwise, we subdivide C into smaller domains and try again. Eventu- 
ally, we succeed in partitioning the problem into six-dimensional pieces, which were 
verified by interval methods. 

If the quad cluster is a mixed case, then [F] gives 

a < vor , — l.OA pt, 

so also 

3 1 
a < -vor +-(-1.04pf). 

In the pure Voronoi case with no quarters and no enclosed vertices, we have the 
approximation 

a < vor(-, \/2) < voro . 

If we prove vor < a(dihi +dih 2 ) + b, the mixed case is established. This is how 
Inequality 4.2.1 was established. Dimension reduction reduces this to a seven- 
dimensional verification. We draw the shorter of the two diagonals between corners 
of the quad cluster. As we begin to subdivide this seven-dimensional domain, we 
are able to separate the quad cluster into two simpliccs along the diagonal, each 
scored by voro. This reduces the dimension further, to make it accessible. The two 
cases, 4.2.2 and 4.2.3, are similar, but we establish the inequalities 

3 1 

- vor + - (- 1 .04 pt) < a(dihi + dih 2 ) + b, 
vor(-, V2) < a(dihi +dih 2 ) + b. 



This completes our sketch of how the verifications were made. 



